Doping induced metal- insulator phase transition in NiO 
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The insulator to metal phase transition in NiO is studied within the framework of reduced density 
matrix functional theory and density functional theory. We find that the spectral density obtained 
using reduced density matrix functional theory is in good agreement with experiments both undoped 
as well as doped NiO. We find that the physical description of the hole-doping induced phase 
transition qualitatively differs depending on whether NiO is calculated within density functional 
theory or reduced density matrix functional. In the former case the underlying mechanism of the 
phase transition is identified to be a rigid shift of chemical potential, while in the latter case a 
redistribution of the spectral weight drives the transition. These latter results are found to be in 
good agreement with both experiments and previous many-body calculations. 

PACS numbers: 



The antiferromagnetic Mott insulators present one of 
the outstanding challenges to the first principles treat- 
ment of extended solids. The prototypical example of 
such a system is NiO, an antiferromagnetic Mott insula- 
tor with a measured gap of 4.1eV and a magnetic mo- 
ment of 1.7 ^LB- The standard local density approxima- 
tion (LDA) within density functional theory (DFT) pre- 
dicts a metallic spectrum for NiO, in fundamental dis- 
agreement with experimental reality. The inclusion of 
spin polarization via the local spin density approxima- 
tion (LSDA) results in a small gap, and a description of 
NiO as a Slater insulator. However, both the gap and 
the magnetic moment are severely underestimated sug- 
gesting that the Slater antiferromagnetic state obtained 
within the LSDA does not describe the true nature of 
NiO. 

In order to over come this deficiency of the LSDA 
Kohn-Sham (KS) spectra, Rodl et al. proposed the use 
of two separate fitting parameters: an on-site Coulomb 
term U and a scissors shift A, where A is the difference 
between the experimental gap and the KS gap obtained 
using LSDA+[/ functional [1]. With a certain choice of 
these two parameters the KS spectra of of NiO can be 
made to agree with computationally expensive many- 
body techniques such as dynamical mean field theory 
(DMFT) PHH] , reduced density matrix functional theory 
(RDMFT)[6], and the GqWo method[l^. This scissors 
corrected LSDA-t-[/ method comes under the heading of 
the so called correlated band theory method. 

What makes NiO even more interesting is its behavior 
as a function of doping: one finds a insulator metal phase 
transition (IMT) on doping the system with Li, which 
amounts to hole doping [7]. The physics of this phase 
transition found to be very rich due to a subtle interplay 
of charge transfer and Mott localization: despite being a 
text book Mott insulator, NiO also has a strong charge 
transfer character due to the large overlap (in energy) 
of the Ni d and 0-p states. Any theory attempting to 
capture this IMT in NiO must be capable of treating 



Mott correlations and charge transfer effects on an equal 
footing, presenting a significant theoretical challenge. 

In the present work we study the IMT in NiO using 
three different approaches: the local spin density ap- 
proximation (LSDA) within DFT, correlated band the- 
ory method, LSDA-I-C/, and a many-body technique, 
RDMFT. In doing so we demonstrate that even though 
at zero doping all these methods give similar spectra, the 
physics of hole doping phase transition is qualitatively dif- 
ferent for the different theoretical methods: within band 
and correlated band theory methods the metalization oc- 
curs due to rigid shifting of the chemical potential into 
the valence band, with the separation between the Hub- 
bard bands remaining approximately constant. In total 
contrast to this, within RDFMT one finds that the phase 
transition is driven by a transfer in spectral weight from 
the upper and lower Hubbard bands to a low energy peak, 
known as the correlated peak. These latter results we 
find to be in excellent agreement with previous many- 
body results obtained using DMFT[S]. 

Recently, RDMFT has shown potential for correctly 
treating band as well as Mott insulators [HI [HI or in 
other words, treating both Mott correlations and charge 
transfer effects on an equal footing. Within RDMFT, the 
one-body reduced density matrix (1-RDM) is the basic 
variable [TOl E] 
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where ^' denotes the many-body wave function and N is 
the total number of electrons. Diagonalization of 7 pro- 
duces a set of orthonormal Bloch functions, the so called 
natural orbitals[TD]. Lpi\^, and occupation numbers, n^k- 
Extending RDMFT to the truly non-collinear magnetic 
casejGi], by treating the natural orbitals as two component 
Pauli-spinors, leads to the spectral representation 
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The necessary and sufficient conditions for ensemble A^- 
representability of 7 were provided, in a classic work, by 
Coleman|12j. These conditions require 
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This method is known to produce accurate spectra for fi- 
nite systems[T5] as well as solids|5]. Following the above 
procedure the spectral density for doped and undoped 
NiO is calculated using the full-potential linearized aug- 
mented plane wave [IB] code Elk|17j. 



In terms of 7, the total ground-state energy [TT] of the 
interacting system is (atomic units are used throughout) 



^[7] = - ^tr, j Jim^ Vl^{v, r') + j p(r)K=xt(r) d 



r - r' 
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where p{r) — trCT7(r, r), Vcxt is a given external poten- 
tial, and i?xc we call the exchange-correlation (xc) energy 
functional. In principle, Gilbert's [TT] generalization of 
the Hohenberg-Kohn theorem to the 1-RDM guarantees 
the existence of a functional E['y\ whose minimum, for 
fixed a T4xt, yields the exact 7 and the exact ground-state 
energy. In practice, however, the correlation energy is an 
unknown functional of 7 and needs to be approximated. 
While there are several known approximations for the xc 
energy functional, the most promising for extended sys- 
tems is the power functional[8| where the xc energy is 
given by 
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with a indicating the power in the operator sense. In 
view of the universality of the functional i?xc[7], the 
value of a should, in principle, be system-independent. 
A few " optimum values" of a have been suggested in the 
literature [H [13 [H]; in the present work a is fixed to 
0.656. 

In order to the study the doping dependent metaliza- 
tion of NiO one crucially needs spectral information. To 
obtain this information from RDMFT, which by its very 
nature is a ground-state theory, is a difficult task. In 
this work we extract spectral density from RDMFT us- 
ing the method recently proposed in Ref. [6]. Within 
this method the diagonal of the spectral density, the so 
called density of states, is determined using the following 
relation: 



DOS = 2ttY, hk<5(w - e'k) + (1 - "*k)<5(w + 4)] (6) 
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FIG. 1: (Color online) Projected and total density of states 
for NiO as a function of energy (in eV). Results are shown for 
undoped (lower panels) and hole doped (upper panels) NiO. 
The results are obtained using two different functionals within 
DFT; the LSDA (left panels) and the LSDA+U method (right 
panels). Experimental data from Ref[TS]is also presented for 
comparison. 



We first examine the behaviour of undoped NiO within 
band and correlated band theory by using LSDA and 
LSDA-|-[/ methods; results for total DOS and the site 
and symmetry projected spectral density are presented in 
Fig. 0. One can observe that while the LSDA band gap 
is grossly underestimated as compared to experiment, a 
rigid shift of the valence band to lower energies would 
result in a spectrum with an overall shape in good agree- 
ment with experiment. Indeed, the correct ordering of 
the t2g and Cg states is obtained within LSDA, i.e., the 
band gap, even though underestimated, separates a 
valence and a eg conduction band, with a substantial 
overlap of eg and t2g states in the valence band. 

On applying an on-site Coulomb repulsion the the gap 
opens, but at the cost of a deterioration in the agree- 
ment of the overall spectral weight with experiment. In 
particular, the band ordering is now found to be incor- 
rect at precisely the value of U that leads to the correct 
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Undoped Doped 



LSDA 




5.42 


4.81 




2.43 


2.43 


O-p 


1.73 


1.51 


LSDA+f/ 




5.52 


5.52 




2.19 


1.80 


O-p 


1.78 


1.51 


RDMFT 




5.49 


5.49 




2.27 


2.27 


O-p 


1.88 


1.08 



TABLE I: Number of electrons in various site and symme- 
try projected states. Results obtained for undoped and hole 
doped (with 1.2 holes per formula unit) NiO obtained using 
the LSDA, the LSDA-|-[/ method, and RDMFT. 



magnitude of the gap, which occurs between purely eg 
states. This striking defect of the LSDA-f [/ treatment of 
NiO was also noticed in Ref. jl, and the use of a smaller 
U with an additional external parameter, the so called 
scissors correction A, was suggested as a remedy. How- 
ever, while the spectrum of the equilibrium ground state 
is improved by this procedure, albeit with the use of an 
additional fitting parameter, such an approach cannot be 
used to study the insulator to metal transition in NiO: 
the very meaning of the parameter A is lost once the 
material enters the metallic phase. 

Upon doping the KS spectrum shows a rather simple 
behaviour for both LSDA and LSDA-}-[/ methods; the 
chemical potential rigidly moves into the valence band 
leading in consequence to mctalization. We may esti- 
mate the distribution of the added hole by integrating 
the site and symmetry projected DOS up to the chemical 
potential. The values of the electronic charge in various 
symmetry states are presented in Table I. It is clear that 
in total contrast to the experimental situation ^ , within 
the LSDA and LSDA+C/ methods the added holes are 
equally distributed between the Ni d and O p-states. 

We now consider the situation where NiO is treated 
within RDMFT, see Fig. Firstly, it is obvious that 
the total spectral density obtained for undoped NiO is 
in good agreement with experiments jl8j. The spectrum 
at high energies (12-15eV), which within DMFT is some- 
what broadened due to analytic continuation, is brought 
out very clearly by RDMFT. The value of the band- 
gap obtained using RDMFT (4.88eV) is larger than ex- 
periments (4.1eV), and the magnetic moment (1. SB/is) 
smaller than the experimental value of 1.7 fj,B- There are 
two reasons for the smaller value of the magnetic moment 
within RDMFT as compared to experiment. Firstly, the 
calculations are performed with the FP-LAPW method 



in which space is divided into spheres around the atoms, 
the so called muffin-tins, and an interstitial region. In 
the case of fully non-coUinear magnetic calculations the 
magnetic moment per site is calculated by integrating 
the magnetization vector field inside the mufhn-tin. This 
implies the loss of a small part of the moment to the in- 
terstitial region. Secondly, the power functional induces 
a slight non-coUinearity in the magnetization leading to 
yet more loss in the integrated z-projected moment. 

The site and symmetry projected spectral density is 
also shown in Fig. (|2| . In good agreement with previous 
many-body studies ^1 [20], we find that the conduction 
band is purely eg in character. A substantial overlap of 
Ni d and Oxygen p states may be seen in the valence 
spectrum, highlighting the presence of charge transfer 
effects in NiO. 

Turning to the hole doping of NiO, we find that the ef- 
fect on the spectral density is strikingly different within 
RDMFT, as compared to both the LSDA and LSBA+U 
methods. Hole doping is found to lead to a redistribution 
of spectral weight from the upper and the lower Hubbard 
bands towards the chemical potential, which in turn leads 
to an IMT. If one uses the correlated band theory defini- 
tion of U as being roughly equal to the distance between 
the upper Hubbard band and the correlated peak, then 
it is evident that the value of U changes as a function of 
dopiirg. An additional striking difference between DFT 
and RDMFT results is that the added hole is mostly lo- 
calized in the O p-states (the precise amount of charge in 
various orbitals is presented in Table I), in marked con- 
trast to the equal distribution between Ni d- and O p- 
states found within LSDA and LSDA-I-C/ method. These 
results are in excellent agreement with experiments in 
which it is observed that the added holes have mainly 
O-p character. 

As shown in Fig. [2j the spectral density of doped 
NiO is also in good agreement with available experimen- 
tal data[7j. Turning to the quantitative description of 
the phase transition afforded by RDMFT, we find that 
the substantial metalization occurs at a somewhat higher 
value of hole doping than observed in experiments (0.5 
holes). There are two principle reasons for this. Firstly, 
the undoped gap for NiO (4.88eV) is larger than the ex- 
perimental value, and hence additional hole doping with 
be required to drive the material to the metallic state. 
Secondly, we do not study the effect of an actual impurity 
added to the system but rather the hole doping is simu- 
lated, as in previous works, by the removal of electronic 
charge from the unit cell while adding a constant com- 
pensating positive background to ensure charge neutral- 
ity. This method is commonly known as the virtual crys- 
tal approximation. Interestingly, previous DMFT calcu- 
lations also identified the transfer in spectral weight as 
the driving mechanism for the hole doping transition in 
NiO[5], and many of the features identified in this arti- 
cle were also found in the DMFT treatment of NiO. We 
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thus find that RMDFT treatment of NiO is much close namical mean field theory, 
to many-body theories such as the DMFT, than to the 

LSDA or correlated band approaches. 
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FIG. 2: (Color online) Total and projected spectral density for 
NiO, obtained using RDMFT, as a function of energy (in eV). 
Results are presented for undoped (lower panel), hole doped 
with 1.2 holes per formula unit (upper panel). Experimental 
results from Ref. 1181 and [71 are presented for comparison. 



To summarize we have presented the physics behind 
the doping induced insulator to metal phase transition in 
NiO and have found that the physics of phase transition is 
bought out in strikingly different ways by different theo- 
retical methods: within DFT based studies metalization 
occurs due to a rigid shifting of the chemical potential 
into the valence band, with the separation between the 
Hubbard bands remaining approximately constant. In 
contrast, within RDFMT the phase transition is driven 
by a transfer in spectral weight from the upper and lower 
Hubbard bands to what is known as the correlated peak. 
The RDMFT results closely resemble those of the dy- 
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